Effective Hamiltonian approach to adiabatic approximation in open systems 
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The adiabatic approximation in open systems is formulated through the effective Hamiltonian 
approach. By introducing an ancilla, we embed the open system dynamics into a non-Hermitian 
quantum dynamics of a composite system, the adiabatic evolution of the open system is then defined 
as the adiabatic dynamics of the composite system. Validity and invalidity conditions for this ap- 
proximation are established and discussed. A High-order adiabatic approximation for open systems 
is introduced. As an example, the adiabatic condition for an open spin-| particle in time-dependent 
magnetic fields is analyzed. 
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I. INTRODUCTION 

As one of the oldest theorems in quantum mechan- 
ics, the adiabatic theorem^] tells us that if a state is an 
instantaneous eigenstate of a sufficiently slowly varying 
Hamiltonian H{t) at one time, then it will remain close to 
that eigenstate up to a phase factor at later times, while 
its eigenvalue evolves continuously. The adiabatic theo- 
rem underlies the adiabatic approximation scheme, and 
has potential applications in several areas of physics such 
as the Landau-Zener transition in molecular physicsQ, 
quantum field theoryQ and geometric phaseQ. Re- 
cently, there has been a growing interest in the adiabatic 
approximation in the context of quantum information, 
for example, geometric quantum computation!^ IS 13 
and the new quantum algorithms 0,113 based on the adi- 
abatic approximation. 

The concept of adiabaticity has been put forward ^ 
and applied to quantum information^^! by Sarandy and 
Lidar in open systems. In their approach, the open sys- 
tem is described by a master equation ip{t) = C{t)p{t) 
= [H{t),p{t)] + L{t)p{t), and the adiabatic approxima- 
tion is characterized by independently evolving Jordan 
blocks, into which the dynamical superoperator C{t) can 
be decomposed. This extension for the adiabaticity in 
open systems is in a systematic manner, though it is chal- 
lenging to calculate the superoperator and determine the 
Jordan decomposition in practice, especially for compli- 
cated disturbances L{t)p{t). The other extension for the 
adiabaticity has been presented by Thunstrom, Aberg 
and Sjoqvist for weakly open systems "13]. In this ap- 
proach, the eigenspace of the system Hamiltonian are 
primary instead of instantaneous Jordan blocks. Thus 
the adiabatic approximation is irrespective of the form 
of L{t)p{t). The main difference between the two ap- 
proaches is that Sarandy et al. focus on the eigenspace 
of the entire superoperator C(t), while Thunstrom et al. 
emphasize the decoupling among the eigenspaces of the 
system Hamiltonian H{t). 

In this paper, we present a new approach to the adia- 



batic approximation in open systems by introducing an 
ancilla to couple to the quantum system. An effective 
Hamiltonian which governs the dynamics of the compos- 
ite system (quantum system plus ancilla) is derived from 
the master equation. We then define the adiabatic limit 
of the open system as the regime in which the composite 
system evolves adiabatically. The validity and invalidity 
conditions for this approximation are also given and dis- 
cussed. An extension to high-order adiabatic approxima- 
tion is presented. As an example, the adiabatic evolution 
of a dissipative spin-i particle driven by time-dependent 
magnetic fields is analyzed. 

The structure of this paper is organized as follows. In 
Sec. II we embed the open system dynamics into a non- 
Hermitian quantum dynamics by introducing an ancilla. 
The adiabatic approximation is defined, and the validity 
condition is given in Sec. HI. In Sec. IV, we present an 
example to get more insight on the adiabatic evolution in 
open systems and discuss the validity condition. Finally 
we present our conclusions in Sec.V. 



II. EFFECTIVE HAMILTONIAN DESCRIPTION 
FOR OPEN SYSTEMS 

We begin with the master equation for the density ma- 
trix p(t)[lj|) 



d 



i-p{t) = [H{t),p{t)]--Y,{Ll{t)Lu{t)p{t) 



+ p{t)Ll{t)Lu{t)-2Lu{t)p{t)Ll{t)} 

= [H{t),p{t)]+L(t)p{t) 

^ mPit), (1) 

where H{t) is a Hermitian Hamiltonian and L(t)p{t) 
describes system-environment couplings and the result- 
ing irreversibility of decoherence, Lk{t) may be time- 
dependent operators describing the system-environment 
interaction. This master equation is of the Lindblad 
form, and thus it guarantees the conservation of Trp(t) 
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and the positivity of probabilities. Our aim in this Sec- 
tion is to solve Eq.(^ by introducing an ancilla cou- 
pling to the system. This approach was first proposed 
in[l5j and called effective Hamiltonian approach. The 
idea of this approach is as the following. The density 
matrix p{t) of the open system can be mapped onto a 
pure state by introducing an ancilla. The dynamics of 
the open system is then described by a Schrodinger-like 
equation with an effective Hamiltonian that can be de- 
rived from the master equation. In this way the solu- 
tion of the master equation can be obtained in terms 
of the evolution of the composite system by converting 
the pure state back to the density matrix. To proceed, 
we assume that the dimension of the system Hamilto- 
nian H{t) is independent of time t, so that we may write 
Hit) = J2n Enit)\Enit)){En{t)\, whcrc TV is a constant 
representing the dimension of the system. The ancilla la- 
belled by A is introduced as same as the open system in 
the sense that its Hilbert space spanned by {|e„i(0))} has 
dimension N and remains unchanged in the dynamics. 
Thus {\En(Q)) (E) \era{0))} may be taken as an orthonor- 
mal and complete basis for the composite system. A 
pure state for the composite system in the iV^-dimension 
Hilbert space may be constructed as 

N 

l^p)- E P^nn{t)\EmmKm, (2) 
m,n— 1 

I 



where Pmn{t) are density matrix elements of 
the open system in the basis {|i?m(0))}, i.e., 
P^n[t) = {E„,{0)W)\Enm. Clearly, p{t)\^ p[t)) = 
Tr{p^{t)) < 1, so this pure bipartite state is generally 
not normalized except that the initial state of the open 
system is pure and the evolution is unitary. With these 
definitions, we now try to find an effective Hamiltonian 
Hrit), such that the bipartite pure state \'^p{t)) satisfies 
the following Schrodinger-like equation 



ijr\^p{t))=nT{mp{t)). (3) 



To simplify the derivation, we write the master equation 
Eq.lUl as 



k 

(4) 

with n{t) = H{t) - i/2J2kLi{t)Lk{t). Substituting 
equation Eq.® together with Eq.|0} into Eq.©, one 
finds, 



z-|*p(t)} = (i^m(0)|HW|i?p(0))(i?p(0)|p(i)|i?„(0))|i?„,(0))|e„(0)} 

- {E^mPit)\Ep{0)){EpiO)\H\t)\E„{0))\E^^iO))\e„{0)) 

+ *E E {En^mLk{t)\Ep{0)){Ep{OMt)\E,{0)){E,{0)\Ll^^^^ 

k m,n,p,q 

= E [WWpW|£;n(0))|e„(0)) -H^(t)pW|i?„(0))|e„(0))] + Lt{t)Lk{t)p{t)\E„{0))K{0)) 

n k,n 

EE nrimpit)), (5) 



where Tixit) is defined by 

HtW - nit) - H^{t) + iY,Lt{t)Lk{t), (6) 

k 

and referred as the effective Hamiltonian. Operators 7i(t) 
and L}^{t) are for the open system, which take the same 
form as in Eq.l^J, while TL^it) and L^{t) are operators 
for the ancilla and defined by 

(e™(0)|6-^|e„(0)) = (ii;„(0)|6^|i?™(0)), (7) 

with O = 7i(t), or Lk{t). The first two terms in the ef- 
fective Hamiltonian Hrit) describe the free evolution of 
the open system and the ancilla, respectively, while the 
third term characterizes couplings between the system 



and the ancilla. The other works aiming at setting the 
passage from pure to mixed states in an unitary evolution 
scheme can be found in Rcf . 0, 0| , but they proceed 
differently. In Ref. J.6J the authors have applied the tech- 
nique of operator splitting to deal with weakly dissipa- 
tive systems. The unitary integrator for the Hamiltonian 
evolution and the conventional integrator for the dissi- 
pation were combined to evolve the open system. Based 
on the wave operator p defined by the density matrix 
p, p = pp\ the unconvention al q uantum mechanical for- 
malism was proposed in Ref.[l7| to study the dynamics 
of open systems. This scheme allows a generalized uni- 
tary evolution between pure and mixed states, and sheds 
new light on the connection between symmetries and con- 
servations laws. Another work[l^ deals with open sys- 
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tems by embedding elements of the density matrix in a 
higher-dimensional Liouville-Bloch equation. The dissi- 
pation and dephasing in the open system were included in 
the non-Hermitian superoperator. Compared with these 
schemes, our effective Hamiltonian approach has the ad- 
vantage that it is easy to calculate, and as will be shown 
in the next section, the extension for adiabaticity from 
closed systems to open systems is straightforward. 



III. THE ADIABATIC APPROXIMATION IN 
OPEN SYSTEMS 

In this section, we introduce an adiabatic approxima- 
tion for open systems. Conditions for this approxima- 
tion are also derived and discussed. We will restrict our 
discussions to systems where the effective Hamiltonian 
TLrii) is diagonalizable with nondegener ate eigenvalues. 
For further details we refer the reader to Ref. [13] , where a 
general discussion on non-Hermitian quantum mechanics 
is presented. 

Let us first define the right and left instantaneous 
eigenstates of TLt (t) by 



nT{t)\R,n{t)) = \,n{t)\Rrnit)), 



(8) 



It is easy to show from Eq.® that {L„i{t)\Rn{t)) ~ (5„m 
for Am 7^ A„. Now we are ready to define the adiabatic 
evolution for open systems, which is directly follow-up 
from that for closed systems. An open system govern 
by the master equation i?g. 0) is said to undergo adi- 
abatic evolution if the composite system govern by the 
effective Hamiltonian TCrit) evolves adiabatically. In 
other words, if the effective Hamiltonian is changed suf- 
ficiently slowly, then the composite system in a given 
non-degenerate eigenstate of the initial effective Hamil- 
tonian TiriO) evolves into the corresponding eigenstate 
of the instantaneous Hamiltonian Ti.T{t), without making 
any population transitions. This leads to the definition 
of adiabatic evolution in the open system which is the 
corresponding adiabatic dynamics in the composite sys- 
tem. This definition is a straightforward extension of the 
idea of adiabatic evolution for open systems, and it will 
be shown below that the condition for this to occur backs 
to the conventional adiabatic evolution when the system 
is a closed system. Let us now derive the validity condi- 
tions for open system adiabatic evolution. To this end, 
we expand |^p(i)} for an time t in the instantaneous right 
eigenstates of Ht (t) as 



I*pW) 



.We" 



''I>-^^^'^\R„,{t)). (9) 



Substituting in into the Schrodingcr-like equation 

Eq.Q, one obtains 

C™(t) + {L^it)\Rrn{t))Crn{t) = Coff{t), (10) 









T 


Time ^ 


P„(0 

T ime ^ 







FIG. 1: In closed systems, the transition rate among distinct 
instantaneous eigenstates {|i5m(t))}of a varying Hamiltonian 
is zero in adiabatic evolution(left), whereas it is zero among 
distinct density matrices {pm(f)} in open systems (right). 
The distinct density matrices pm{t) can be obtained by map- 
ping the instantaneous eigenstates of the effective Hamilto- 
nian back to the open system. 



where 

Coffit) c4t)(imWl^n(t))e-'/o(^"(^)-^"(^»'^^ 

(11) 



Formal integration of Ea. HlU|) yields. 



c„(t) = Cm(0) 



{LmiT)\RmiT))Cm{T)dT 



Coff{T)dT. 



(12) 



In accordance with the definition of adiabaticity in open 
systems, the adiabatic regime is obtained when Co//(t) = 
Cof f{T)dT is negligible. This condition ensures that 
the mixing of coefficients Cm (t) corresponding to distinct 
eigenvalues Xm{t) is absent, which in turn guarantees 
that the change in T-Crit) is sufficient slow. The latter 
claim can be shown by rewriting the adiabatic condition 
Coffit) ^ as 



{Ln{t)\dHT{t)/dt\R^{t)) 



(A„(t)-A„(t))2 
(L„(t)|7?™(t)) 



Xmit)-Xn{t) 



< 1. 



(13) 



In terms of the density matrix p{t), (i„(t)|i?„i(i)) 
can be expressed as {Ln{t)\Rm{t)) = Tr{pn{t)pm{t)), 
where the elements of Pm(i) are defined as [pmit)]a0 = 
{ei3{0)\{Ea{0)\Rm{t)) , i.e., Pm{t) is the density matrix 
corresponding to the m-th right eigenstate of the effec- 
tive Hamiltonian. Then the adiabatic condition in this 
case becomes \Tr{pn{t)pm{t})\ ^ \Xm{t) — A„(t)|. This 
relation implies that the transition rate from one path 
Pn{t) to the other path Pm{t) is negligible in the adi- 
abatic evolution. From the other aspect, the adiabatic 
evolution for open systems indicates that \{Lm{t)\'9 p{t))\ 
remains constant in the dynamics for any m (degenerate 
levels are excluded). This results in Tr{pm{t)p{t)) = con- 
stant, showing again that transition rate among distinct 
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paths Pm{t) is zero in the adiabatic dynamics. A compar- 
ison of the adiabatic evolution in closed systems with that 
in open systems was sketched in figure ^ For a closed 
system, L{t)p{t) = 0, yielding Ht(0 = H{t) - H^{t), 
where H{t) is the system Hamiltonian given in Eq.(^, 
H^{t) is the counterpart of H{t) for the ancilla. Note 
that H^{t) takes the same form as H{t), except that it 
operates on the ancilla. Then the adiabatic condition in 
this case backs to 



{E^{t)\E^{t)) 



< 1. 



(14) 



This is the well known adiabatic condition for closed 
system[23. In other words, the adiabatic condition de- 
fined here for open systems back to the adiabatic condi- 
tion for closed systems when L{t)p{t) = 0, showing the 
consistency of the definition for open systems. 

The definition of adiabatic evolution can be easily gen- 
eralized to high-order adiabatic approximation. To this 
end, we take the instantaneous eigenspace of the system 
Hamiltonian H{t) as the primary Hilbert space. As will 
be shown, this choice has the advantage that the for- 
malism would return straightforwardly to Thunstrom's 
results for weakly open systems. Defining 



where H{t) = H{t) + Z{t) - ^J2k Ll{t)Lk{t). The high- 
order adiabatic evolution is then defined as the adia- 
batic evolution of the composite system with the effective 
Hamiltonian Ti.T{t). The reason of referring the adiabatic 
approximation with Tirit) as the hight-order adiabatic 
approximation is that Z{t) contains off-diagonal terms 
Em^n \Em{0)){E^{t)\En{t)){En{0)l m + n. In the adi- 
abatic approximation defined with HTit), these terms 
have been ignored when the open system approaches to 
a closed system. However, they could not be ignored in 
the latter definition. In fact, the adiabatic approximation 
defined with Hrit) is for p{t) in the master equation H17|l . 
which is different from p{t) in Eq.(^, and this definition 
is of consistency with the high-order approximation in 
closed systems. This point can be found by the same 
analysis presented above for the definition with Tirit)- 



IV. EXAMPLE: THE ADIABATIC EVOLUTION 
OF AN OPEN SPIN-i PARTICLE IN 
TIME-DEPENDENT MAGNETIC FIELDS 



N 



t/(t)=5]|i?„(0))(i?„(i)|, 



(15) 



and 



pit) = U{t)pit)UHt), (16) 
one gets a master equation for p{t) from Eq.lQJ I^Ij 



i-m = [Hit),p{t)] + [zit),m] 



Y,{Li(t)h{t)m + p{t)Ll{t)Lk(t) 



where 



2U{t)p{t)Ll{t)}, 



H{t) = Y,En{t)\E,,{0)){Enm 



(17) 



In this section, we present an example to get more 
insight of the adiabatic evolution in open systems. The 
example consists of a dissipative spin-i particle driven by 
a time-dependent magnetic field. The master equation 
govern the dynamics of such a system can be written as 



-i[H{t), p] + -{2a-_po-+ - pcr+cT- - cr+cr^p}, (19) 



where H{t) — pB{t) ■ a denotes the system Hamilto- 
nian, K stands for the spontaneous emission rate. (Ti,i — 
z,+,— are Pauli matrices, and CTz = |e)(e| — \g){g\ (|e), 
denotes the state of spin-up, and \g) spin-down). Sup- 
pose B{t) = So (sin 6* cos 0, sin 61 sin 0, cos 6*), the effective 
Hamiltonian corresponding to p{t) reads. 



Z{t) = iU{t)uHt), Lkit) = U{t)Lk{t)U\t). 

The effective Hamiltonian corresponding to this master 
equation can be derived by the same procedure and given 
by 

nT{t) = nit) -n^{t) + iJ2 Lk{t)Lk{t), (is) 

k 



HT{t)=n{t)~'H^{t)+iK(T^T^, (20) 



where n{t) = H{t) - i|e)(e|, H^{t) = H^{t) + ||e)A(e|, 
and T_ — \g)A{e\ represents the Pauli matrix of the an- 
cilla. In a subspace spanned by {\gg), \ge), \eg), |ee)}, the 
effective Hamiltonian can be written as (in units of pBo), 



5 



sin 6*6*'^ 

— Z0 



sm f e 




'2cose - 0.5^7 


sin 6*6"'"^ 



sin 



2 008 6* - 0.5^7 

-sinSe'"^ 



Z7 



smf 



(21) 



-27 



where 7 = k/ fiBo, and the first (second) letter in the basis 
\ab){a,b = 6,(7) stands for states of the system(ancilla). 
The eigenvalues Xj (in units of ^Bq) of Hrit) are given 
by 

Ai = 0, 

{Xj + O.bijf + 0.5ijiXj + 0.5^7)2 - 4(Aj + 0.5^7) 



= 2z72cos2 0, (j^ 2,3,4), 
and the corresponding right eigenstates 



1^.) = 



1 



as well as the left eigenstates 
1 



Here 



'-3 
ai 



B, 



(22) 



(23) 



(24) 



4 cos^ 6* + 0.2572 



= 1 



sin^ 6* 
2 cos 61 - 0.57^6'"^ 
sin 6* 

2 cos 61 + 0.57^6-*'^ 



hi = 

ci = 

smw 

di = 1; 

Ai = £>! = = Ci 0. 



For J = 2,3,4, 



(25) 



A, 



-(2 cos 6* + 0.5i7 + A^), 
2 sin 6le''^, 

2 sin 61(2 cos 6* + 0.5i7 + A-,)e 
2 cos6' — 0.5z7 — Aj 

~{2cos0-O.5i-f + Xj), 
A, 



-i(p 



h - Aj 

i7 + A 



(Dj-Aj) 



2 cos 6* + 0.5^7 + Xj ' 

sin6ie*'^(Dj ~ Aj) 
2 cos 9 — Q.5ij — Xj ' 




FIG. 2: (Color online) An illustration of F = F(7,aj) as a 
function of uj and 7. lj was chosen in units of fiBo, and this 
plot is for 9 = jTT. 



To simplify the discussion, we assume that cj) — tot with 
a constant w, and remains unchanged in the dynam- 
ics. With this assumption, it is readily to show that 
all eigenvalues of Hrit) are time-independent, so Mj are 
constant. To show the dependence of the adiabatic condi- 
tion on the spontaneous emission rate 7 and w, we define 
the following function with max taken over all m and 



r(7, uj) = max ■ 



{Ln{t)\Rra{t)) 



Xr 



An 



1,2,3,4, 



(26) 



(27) 

which characterizes the violation of the adiabatic evolu- 
tion. The numerical results of r(7,<x') versus 7 and uj 
were illustrated in figure 1. The rotating frequency uj of 
the driven field was plotted in units of ^Bq, where Bq 
is the modulus of the driving field. Figure 1 shows us 
that r(7,[jj) depends on uj linearly with a fixed 7, while 
r(7,tt') decreases with 7 increasing. This can be under- 
stood as the following. Without the driving field B{t), 
the dependance of \Rm{t)) on 7 behaves like e~^*, while 
Am depends on 7 linearly. Therefore, r(7, w) decays with 
7 increasing and tends to zero when 7 — > c». 



V. SUMMARY AND DISCUSSION 

We have presented an adiabatic approximation scheme 
for open systems. In contrast with the conventional adi- 
abatic approximation, the adiabatic approximation for 
open systems have been defined as the adiabatic dynam- 
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ics of a composite system, which consists of the system 
and an ancilla. Our effective-Hamiltonian based defini- 
tion of adiabaticity retains the conventional adiabatic ap- 
proximation in the ideal case of closed systems, hence it 
is of consistency. The definition of adiabaticity in open 
systems has been extended to a high-order adiabatic ap- 
proximation, which was defined in accordance with an 
master equation for the rotated density matrix. The 
validity and invalidity condition of adiabticity has been 
derived and discussed. The violation of adiabatic evolu- 



tion has been demonstrated by an field-driven dissipative 
spin-i particle. The other applications and demonstra- 
tions such as geometric phase in open systems and the 
effect of decoherence on quantum adiabatic computing 
will be addressed elsewhere. 
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